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Abstract 

We revise existing limits on the D-dimensional Planck scale Mo from the nonobservation of 
microscopic black holes produced by high energy cosmic neutrinos in scenarios with D = 4+n 
large extra dimensions. Previous studies have neglected the energy radiated in gravitational waves 
by the multipole moments of the incoming shock waves. We include the effects of energy loss, as 
well as form factors for black hole production and recent null results from cosmic ray detectors. 
For n > 5, we obtain Mjj > 1.0 — 1.4 TeV. These bounds are among the most stringent and 
conservative to date. 
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I. INTRODUCTION 



Forthcoming colliders fil, cosmic ray observatories 0,0,0], neutrino telescopes 0, and 
space-based experiments jy] will be able to observe black holes (BHs) if the fundamental 
scale of gravity is sufficiently close to 1 TeV Observations of highly characteristic BH 
events at any of these facilities could conceivably provide the first evidence for the existence 
of extra dimensions and make possible the direct study of strong quantum gravity effects 
and strings. On the other hand, the lack of such events in any experiments to date leads to 
lower limits on the scale of higher- dimensional gravity 

To make useful predictions about higher-dimensional gravity based on observations of 
such events, or their absence, a quantitative understanding of the process of BH production 
in high-energy collisions is required. An intuitive picture of this process is provided by a 
simple model known as Thome's hoop conjecture [9(, according to which a BH forms in a 
two-particle collision when and only when the impact parameter is smaller than the radius 
r s of a Schwarzschild BH of mass equal to the total center-of-mass energy -Ecm- The hoop 
conjecture thus predicts a total cross section for BH production equal to the area subtended 
by a "hoop" of radius r s : 

°& = ^ S {E CM ) . (1) 

Up to now, all studies of BH production in TeV-scale gravity have been based on this rather 
heuristic cross section, and have thus been subject to substantial theoretical uncertainties. 

Relatively recently, significant progress has been made in determining the cross section for 
BH production. Early analytic calculations in four dimensions rJ n|| for head-on collisions 



illustrated the process of horizon formation and found that the mass of the final BH was 
about 84% of the initial center-of-mass energy. These calculations were extended to nonzero 
impact parameter by Eardley and Giddings |l2j|, who analytically derived a lower bound on 
the total cross section of approximately 65% of Eq. (0). Relatively recently, a calculation of 
the cross section in higher dimensions was performed by Yoshino and Nambu using numerical 
techniques [r| . In addition, these authors observed significant reductions in the mass of the 
final-state black hole as a function both of impact parameter and dimension. 

If TeV-scale gravity is realized in nature, then the first observational evidence for it will 
likely come from BH-mediated neutrino cosmic ray showers. Ultra-high energy cosmic rays 
hit the Earth with collision center-of-mass energies ranging up to roughly 10 5 GeV. QCD 
cross sections dominate over the BH production cross section by a factor of roughly 10 9 . 
Thus, black holes produced by hadronic cosmic rays are effectively unobservable. This is not 
the case for incoming neutrinos, whose cross section for producing black holes can be orders of 
magnitude larger than SM cross sections, but much less than hadronic [2] . As a consequence, 
neutrinos interact with roughly equal probability at any point in the atmosphere, and the 
light descendants of the black hole may initiate quasi-horizontal showers in the volume of air 
immediately above the detector. Because of these considerations the atmosphere provides 
a buffer against contamination by mismeasured hadrons, allowing a good characterization 
of BH-induced showers when the BH entropy S ^> 1 3]. Additionally, neutrinos that 
traverse the atmosphere unscathed may produce black holes through interactions in the ice 
or water 0. Because the BH production cross section is suppressed by a power of the 
fundamental Planck scale (approaching M D for large numbers of extra dimensions), the 
absence of neutrino showers mediated by black holes implies lower bounds on Mp. 

In this article, we bring up to date existing limits on Mo from the nonobservation of 

BHs at cosmic neutrino detection experiments. Besides incorporating the cross section and 
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energy loss results of Yoshino and Nambu, we also make use of updated parton distribution 
functions and recently available cosmic ray data. 



II. ENERGY LOSS IN BLACK HOLE CREATION 



Previous calculations of the cross section for producing a BH have neglected energy loss 
in the creation of a BH, assuming that the mass of the created black hole Mbh was identical 
to the incoming parton center-of-mass energy \^S. However, recent work |l3[ has shown that 
the energy lost to gravitational radiation is not negligible, and in fact is large for larger n 
and for large impact parameters. The trapped mass (called Ma.h. in Ref. |l3j, and which 
we continue to call M B h [HI), is given by 

M BU (z) = y(z)V§ , (2) 
where the inelasticity y is a function of z = b/b max . Here b is the impact parameter and 



&max = V F(n) r s (V s, n, M D ) (3) 
is the maximum impact parameter for collapse, where 

-2n 7r (n-3)/2p^n+3 
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(4) 



is the radius of a Schwarzschild BH in (A+n) dimensions |16(, and F(n) is the form factor 
explicitly given in Ref. 

This complicates the parton model calculation, since the production of a BH of mass 
Mbh requires that s be M^ R /y 2 (z), thus requi ring the lower cutoff on parton momentum 
fraction to be a function of impact parameter [12(. In what follows we take the uN cross 
section as an impact parameter-weighted average over parton cross sections, with the lower 
parton fractional momentum cutoff determined by the requirement M^ 1 = x^-^M^ [l8| . 
This gives a lower bound x^ lin M| ) /[y 2 (z) s] on the parton momentum fraction x. With this 
in mind, the vN — > BH cross section is 

a(E u ,x min ,n,M D ) = f 2z dz f )2 dx F(n) irrl(V§,n, M D ) fi(x,Q) , (5) 

where x m i n is determined by the requirement that the BH have at least an approximate 
semi-classical description, s = 2xm^E v , i labels parton species, and the fi(x,Q) are parton 
distribution functions (pdfs) 0]. 

The choice of the momentum transfer Q is governed by considering the time or distance 
scale probed by the interaction. Roughly speaking, the formation of a well-defined horizon 
occurs when the colliding particles are at a distance ~ r s apart. This has led to the advocacy 
of the choice Q ~ r" 1 |^|, which has the advantage of a sensible limit at very high energies. 
However, the dual resonance picture of string theory would suggest a choice Q ~ V^- 
Fortunately, as noted in Refs. 0], the BH production cross section is largely insensitive 
to the details of the choice of Q. In what follows we use the CTEQ6M pdfs with 
Q = min{r7 1 , 10 TeV}. 
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FIG. 1: Cross sections a{vN — > BH) for n = 1, . . . , 7 from below, assuming Mb = 1 TeV and 
£min = 1- Energy loss has been included according to Eq. (|SJ). The SM cross section a(vN — > ^X) 
is indicated by the dotted line. 

In Fig. Q we show the BH production cross section for n = 1, . . . , 7 extra dimensions with 
energy loss incorporated as given in Eq. (J3J). The rapid rise in cross section is pushed to 
higher E v than in the case with energy loss neglected. However, the cross sections are still 
well above the SM cross section at E v ~ 10 8 GeV and above, where, as we will see, the 
cosmogenic neutrino flux is large. 

III. COSMIC NEUTRINO DETECTORS 

Energy loss also impacts event rates at cosmic neutrino detectors, not only because the 
cross section is modified, but also because the apertures of cosmic neutrino detectors are 
functions of shower energy. Let Na be Avogadro's number, A(yE u ) the neutrino aperture 
of a given experiment for shower energy yE u , and T be the experiment's running time. The 
number of neutrino showers mediated by BHs is then 

jV(x min ,n,M D ) = N A T J dE u J 2z dz / ( )a dx — A(yE u ) 

xF(n)7rr 2 s (^fs,n,M D )J2fi(XiQ) , (6) 

i 

where d^jdE v is the diffuse flux of cosmic neutrinos hitting the Earth. 

There are several techniques employed in detecting neutrino showers [20]. The most 
commonly used method involves giant arrays of particle counters that sample the lateral 
and temporal density profiles of the muon and electromagnetic components of the shower 
front. Another well-established method involves measurement of the air shower evolution 
- its growth and subsequent attenuation — as it develops by sensing the fluorescence light 
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produced via interactions of the charged particles in the atmosphere. A third method exploits 
naturally occurring large volume Cerenkov radiators such as deep water or ice. Especially 
useful at relevant center-of-mass energies for BH production is emission of Cerenkov radiation 
at radio frequencies. For fluorescence data, a direct measurement of the depth of shower 
maximum X max and the shape of the longitudinal profile provide sensitive diagnostics in 
discrminating between neutrino and hadron showers. In the case of surface arrays, the 
composition information is extracted from a number of shower characteristics which reflect 
the depth of shower maximum and the ratio of muon to electromagnetic content of the 
shower. 

The AGASA Collaboration [21J reports no significant enhancement of deeply-developing 
shower rates given the detector's resolution. Specifically, there is only 1 event observed, 
consistent with the expected background of 1.72 from hadronic cosmic rays. For details, see 
Ref. 0. 

The Fly's Eye detector ceased operation in July 1992 after a life of 11 years. It was 
designed to collect the atmospheric nitrogen fluorescence light produced by air shower par- 
ticles on moonless nights without cloud cover, achieving an overall duty cycle of ~ 10%. 
The experiment recorded more than 5000 events, but no unusual deeply developing showers 
have been found 

Recently, data from an upscaled version of the Fly's Eye experimet have become avail- 
able [2^. The effective aperture of the High Resolution Fly's Eye detector is on average 
about 6 times the Fly's Eye aperture, with a threshold around 10 8 GeV. The instrument 
includes two sites (HiRes I and II) located 12.6 km apart. Each site consists of a large num- 
ber (22 at HiRes I and 42 at HiRes II) of telescope units pointing at different parts of the 
sky. Between November 1999 and September 2001, 1198 events were recorded with at least 
one reconstructed energy greater than 10 8 ' 7 GeV |23j. Because of bad weather conditions, 
272 events were discarded from the sample. None of the 723 events that survived all of the 
cuts required for stereo- mode triggering has A max > 1200 g/cm 2 . Additionally, there are no 
events detected in monocular mode with X max > 1500 g/cm 2 . In the spirit of Ref. (24j |. we 
parametrize the HiRes aperture for deeply (A max > 1500 g/cm 2 ) developing showers by 



A(E V ) = U 



E 

2.7 + log' 



10 10 GeV 



0.5 > km 3 we sr . (7) 



We note that there is an additional small contribution to the HiRes exposure in the energy 
range 10 8 — 10 9 GeV from data collected during 2878 hours livetime [25j. 

The Radio Ice Cerenkov Experiment (RICE) is designed to detect the radio frequency 
Cerenkov radiation produced by neutrino- induced showers in ice j26[. Specifically, the elec- 
tromagnetic channel of the shower produces a radio pulse with a duration of a few nanosec- 
onds and with power concentrated around the Cerenkov angle. Several radio antennae 
positioned in the ice allow for reconstruction of the interaction vertex. For primary energies 
above 10 9 GeV, the Landau Pomeranchuk-Migdal (LPM) effect j2?| leads to a significant 
suppression of the Bethe-Heitler cross sections for the pair production and Bremsstrahlung 
processes in dense materials, and thus dramatically changes the character of the development 
of electromagnetic showers. 

Almost instantaneously after its formation, the TeV-scale BH decays |28j], predominantly 
through radiation of standard model (SM) particles [Hj. About 75% of the BH energy is 
carried off by quarks and gluons and rougly a third of this energy goes into the electro- 
magnetic channel via tt° decay. Only about 5% of particles directly emitted from the BH 
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FIG. 2: The monotonically rising curves are the exposures as functions of shower energy for AGASA 
(solid), Fly's Eye (dotted), HiRes (short dash), and RICE (long dash). The remaining solid curve, 
with a peak around 10 8 ' 5 GeV, is the cosmogenic neutrino flux. 



(z/'s, r's, /x's) do not partake in the shower. The rest of the energy eventually devolves into 
secondary electrom agne tic cascades with particle energies below that for which the LPM 
effect is important [30j. As a conservative estimate we model the aperture for neutrino 
showers mediated by BHs using the hadronic effective volume reported by the RICE Collab- 
oration [jyj with average inelasticity (y) = 0.8. This estimate is supported by the fact that 
BH-induced showers mimic SM neutral current events, characterized by hadronic dominated 
showers with no leading charged lepton. A more rigorous analysis of the BH acceptance at 
the RICE facility is underway |82| . 

The RICE detector comprises 16 dipole radio receivers installed in the holes drilled for the 
AMANDA experiment in the Antaractic ice. Four transmitter antennae are also deployed 
in the ice for calibration purposes. The trigger requires a 4-fold coincidence within a 1.2 /xs 
time window, and various cuts are applied to reject thermal and anthropogenic backgrounds. 
For example, shower vertices are reconstructed using a \ 2 fit to the signal arrival times, and 
the resulting fits are required to be of sufficient quality and to indicate vertices at least 50 m 
below the ice's surface. During a 1 month run in August of 2000 with a livetime of 333.3 
hours, a total of 22 events passed all the automated cuts. These events were then scanned 
for quality and the RICE Collaboration concluded that there are no events consistent with 
neutrino sources |3l| . 

The relative exposures for the different experiments are given in Fig. |21 For details on 
the apertures of AGASA and Fly's Eye, the reader is referred to our previous paper (l4| . 
All in all, there is only 1 event observed with an expected background of 1.72 from hadronic 
cosmic rays, leading to a 95% CL limit of 3.5 BH events 33]. 

To derive the bounds on Mp, we use the "guaranteed" flux of cosmogenic neutrinos arising 
from the decay of ^ produced in collisions of ultra-high energy protons with the cosmic 
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FIG. 3: 95% CL lower limits on the fundamental Planck scale as a function of x m ; n for n = 1, . . . , 7 
extra dimensions (from below). 



microwave background. As in our previous analyses, we conservatively adopt the estimates 
of Protheroe and Johnson [34| with nucleon source spectrum scaling as d&jy/dE oc E~ 2 
and extending up to the cutoff energy 10 12 5 GeV. The total ultra-high energy cosmogenic 
neutrino flux is also shown in Fig. |21 



IV. BOUNDS 



In Fig. El we show 95% CL lower bounds on Mb as derived from Eq. © using the 
exposures and the cosmogenic neutrino flux given in Fig. El requiring M < 3.5 events to be 
observed in cosmic neutrino data samples [35] . The BH entropy is a measure of the validity 
of the semi-classical approximation. For x m i n ^ 3 and n > 5, the entropy 

n + 2 v ' 

yielding small thermal fluctuations in the emission process [36]. Hence, for x m \ n ^ 3 and 
n > 5, strong quantum gravity effects may be safely neglected. Moreover, gravitational 
effects due to brane back-reaction are expected to be insignificant for M B h well beyond 
the brane tension, which is presumably of the order of Mb- The uncertainty illustrated in 
Fig. El associated with x min only concerns BH production and is highly insensitive to decay 
characteristics This is because the signal in neutrino detection experiments relies only 
on the existence of visible decay products. Whatever happens around x m i n ~ 1, it seems 
quite reasonable to expect that BHs or their Planckian progenitors will cascade decay on 
the brane. 

String theory provides a more complete picture of the decay for M B h close to M D , which 
may further justify setting £ min = 1. (Such arguments do not address the issue of brane 
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back-reaction, however.) In string theory, the ultimate fate of the black hole is determined 
by the string/BH correspondence principle [38j]: when the Schwarzschild radius of the black 
hole shrinks to the fundamental string length i s 3> £d, where Id is the fundamental (4 + n)- 
dimensional Planck length, an adiabatic transition occurs to a massive superstring mode. 
Subsequent energy loss continues as thermal radiation at the unchanging Hagedorn temper- 
3^| . The continuity of the cross section at the correspondence point, parametrically in 
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both the energy and the string coupling, provides independent support for this picture [40] . 
Thus, the cross sections given in Fig. [T] can be thought of as lower bounds on o as Mbh 
approaches Mr, [ill ]. 



V. SUMMARY 

Incorporation of the results of Ref. [l3j has eliminated many of the sources of uncertainty 
enumerated in Ref. [H and recapitulated in Ref. [iij. In Ref. |H we identified two sources 
of uncertainty that could reduce the total cross section: the reduction of the mass of the 
final-state BH relative to the initial center-of-mass energy, and expectations for a reduced 
cross section at nonzero center-of-mass angular momentum. 1 On the other hand, we pointed 
out that the classical photon capture cross section and nonrelativistic estimates suggest a 
possible enhancement to the naive geometric cross section nr^ by a factor of 2 or more. (The 
claim of that this upside uncertainty casts doubt on the program of setting limits on 
Md from nonobservation of BHs is mistaken.) Thus we concluded that, in the absence of a 
better quantitative understanding of the process of BH formation, the naive geometric cross 
section provided a reasonable estimate. 

With such calculations now in hand we have repeated our analysis and eliminated 
much of the uncertainty contained in our previous limits on Md, as well as incorporating 
updated exposures from the HiRes and RICE facilities and updated pdfs. (Incidentally, using 
the new pdfs contributed a net difference of about 2% to our results, confirming our previous 
claim 0] that there is very little sensitivity to different choices of pdfs. Furthermore, the 
bulk of the sensitivity is for pdfs at x > m^ H / (2m^E v ) ~ 1CT 2 and large Q, where the pdfs 
are expected to be quite accurate.) In the course of our analysis we observed a competition 
of effects leading to corrections to our previous estimates: enhancement of the geometric 
cross section by form factors of up to 1.9 and enhancement of apertures from new cosmic ray 
data, but a simultaneous reduction in the rate of production of BHs of mass greater than 
x m i n M£>, after taking energy losses into account. It turns out that the latter effect dominates 
and leads to a slight weakening of our limits on Mp. At the same time, our limits are now 
on a much firmer theoretical footing, and maximally conservative in all respects. 

In Fig. E] we compare the bounds derived in this article with existing limits on the fun- 
damental scale of large extra dimensions. Tests of the gravitational inverse-square law at 
length scales well below 1 mm show no evidence for short range Yukawa interactions. For 
n = 2, this negative result can be translated into a 95% CL upper limit of 150 /zm on the 
compactification radius of flat extra dimensions, or equivalently to a unification mass scale 



In fact, a slight modification of our estimate for the modification of the cross section due to angular 
momentum effects has been used to give a surprisingly accurate postdiction of the higher-dimensional 
cross section for BH production 43], providing further evidence for the correctness of the Yoshino and 
Nambu calculations. 
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FIG. 4: Bounds on the fundamental Planck scale Mb from tests of Newton's law on sub-millimeter 
scales, bounds on supernova cooling and neutron star heating, dielectron and diphoton production 
at the Tevatron, and nonobservation of BH production by cosmic neutrinos. The uncertainty in 
the Tevatron bounds corresponds to the range of brane softening parameter 6 (Md/2, Mp); for 
details see Ref. 0. The range in the cosmic ray bounds is for x m i n = 1 — 3. 



M D > 1.8 TeV 0. In such toroidal compactifications the accessibility of towers of Kaluza- 
Klein gravitons may drastically affect the phenomenology of supenovae and neutron stars. 
For n < 3, anomalous cooling of supernovae due to bulk graviton emission and neutron star 
heating by decay of gravitationally trapped Kaluza-Klein modes provide limits on Mb that 
greatly exceed 1 TeV . 

For n > 4 the sensitivity of table-top experiments and astrophysical observations to 
TeV-scale gravity is largely reduced: already for n = 4 (n = 5) supernova cooling yields 
M D > 4.0 TeV (M D > 0.8 TeV) @. For n > 5, the best existing limits on TeV-scale 
gravity are from the absence of trans-Planckian signatures (BH/stringball production) in 
neutrino detection experiments discussed here, and from searches for sub-Planckian signa- 
tures (graviton emission and virtual graviton exchanges) at the Tevatron j4?| and LEP |48| . 
For n > 5 we have derived conservative bounds incorporating the lower limits on the the 
mass trapped in the gravitational collapse. The resulting bounds, M D > 1-0 — 1.4 TeV for 
competitive with those obtained in colliders and among the most stringent 

to date. 
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